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^ ■ Abstract 

■ Fluctuations of charged particle number are studied in the canonical ensemble. 

In the infinite volume limit the fiuctuations in the canonical ensemble are different 
^ , from the fluctuations in the grand canonical one. Thus, the well-known equivalence 

\^ ' of both ensembles for the average quantities does not extend for the fluctuations. In 

lO ' view of a possible relevance of the results for the analysis of fluctuations in nuclear 



o 
o 
o 



X 



collisions at high energies, a role of the limited kinematical acceptance is studied. 



1. Introduction. The statistical approach to strong interactions is surprisingly success- 
ful in describing experimental results on hadron production properties in nuclear collisions 
^ ! at high energies (see e.g. Ref. |1. and references therein). This motivates a rapid devel- 
I opment of statistical models and it raises new questions, previously not addressed in 
statistical physics. In particular, an applicability of the models formulated within vari- 
ous statistical ensembles has been considered for average quantities. The micro-canonical 
ensemble, where the motional and material conservation laws are strictly fulfilled in all 
microscopic states of the system, has to be used for collisions in which a small number of 
particles is produced, like p+p interactions at low energies (see e.g. Ref. j2]). The canoni- 
cal ensemble ( c.e. ), where only the material conservation laws are obeyed, is relevant for 
the systems with a large number of all produced particles, but a small number of carriers 
of conserved charges like electric charge, baryon number, strangeness or charm (see e.g. 
Ref. |3|). Finally, models formulated using grand canonical ensemble ( g.c.e. ) can be 
used when the number of carriers of conserved charge is large enough (see e.g. Ref. ^). 
In the latter approach both material and motional conservation laws are relaxed and the 
mean values of conserved charges and energy are adjusted by introduction of chemical 
potentials and temperature, respectively. 

The question of applicability of various statistical ensembles for the study of fluctu- 
ations of physical quantities has not been addressed up to now. In the text-books of 
statistical mechanics, the particle number fiuctuations are considered in g.c.e only. It 
is because the discussion is limited to the non-relativistic cases, so that in the c.e. the 
particle number is fixed. However, in the relativistic case, relevant for the models of 
hadron production in high energy nuclear collisions, only conserved charges are fixed, and 
consequently the particle number fluctuates in both c.e. and g.c.e. . 
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The analysis of fluctuations is an important tool to study a physical system created 
in high energy nuclear collisions (see e.g. |5j). Recently, rich experimental data on fluc- 
tuations of particle production properties in nuclear collisions at high energies have been 
presented (see e.g. presentations at "Quark Matter 2004"). In particular, intriguing 
results concerning the particle number fluctuations in collisions of small nuclei at the 
CERN SPS have been shown [Bj. These new results motivate our work, in which the 
particle number fluctuations are calculated in c.e. and compared with those obtained in 
g.c.e. . Finally, the possible influence of the limited experimental acceptance on observed 
fluctuations is studied. 

2. Partition function in g.c.e. and c.e. Let us consider the system which consists 
of one sort of positively and negatively charged particles (e.g. vr"*" and vr" mesons) with 
total charge equal to zero (5 = 0. In the case of the Boltzmann ideal gas (the interactions 
and quantum statistics effects are neglected) in the volume V and at temperature T the 
g.c.e. partition function reads: 

z,.uv,T) = J2 J2 -1^-1^ = "^p^^+" + = "^p^2^) • 

N+=ON^=0 +■ 

In Eq. (0) z is a single particle partition function 
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Pdkexp = — Tm'K^[-) , (2) 



where m is a particle mass and K2 is the modified Hankel function. Parameters A+ and 
A_ are auxiliary parameters introduced in order to calculate the mean number and the 
fluctuations of positively and negatively charged particles (the chemical potential equals 
to zero to satisfy the condition {Q) g.c.e. = 0). They are set to one in the final formulas. 

The c.e. partition function is obtained by an explicit introduction of the charge con- 
servation constrain, — = for each microscopic state of the system and it reads: 

ZUV,T) = ± ± ihll}^ S(N, - NJ, = (3) 

N+=ON-=0 +■ 



- / d(P exp [z (A+ e'^ + A_ e"'^)] = Io{2z) . 



In Eq. the integral representations of the 5-Kronecker symbol and the modified Bessel 
function were used: [7j 



5{n) = — / dcj) exp(m0) , Iq{2z) = — / d(j) exp[— + 2z cos0] . (4) 
2vr Jo 27r 



3. Mean particle number. The average number of A^4. and can be calculated as 
and in the g.c.e. ((T)) they are equal to: 



2 



g.c.e. 



(6) 



A± = 1 



By construction the mean total charge is equal to zero: {Q) g.c.e. = g.c.e. — g.c.e. = 
0. In the c.e. Q the charge conservation Q = N^—N_ = is imposed on each microscopic 
state of the system. This condition introduces a correlation between particles which carry 
conserved charges. The average particle numbers are [S]: 
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A± = 1 



Io{2z) 



(7) 



The exact charge conservation leads to the c.e. suppression {Ii{2z) / Iq{2z) < 1) of the 
charged particle multiplicity relative to the result for the g.c.e. ©• The ratio of {N±) 
calculated in the c.e. and g.c.e. is plotted as a function of z in Fig. 1. 
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Figure 1: The ratio of {N±)c.e. © to {N±) 



g.c.e. 



z 

as a function of z. 



In the large volume limit {V ^ oo corresponds also to z ^ oo) the results for mean 
quantities in the c.e. and g.c.e. are equal. This result is referred as an equivalence 
of the canonical and grand canonical ensembles. It can be obtained using an asymptotic 
expansion of the modified Bessel function jTj: 



lim In{2z) 



exp{2z) 



71 Z 



1 - 



16z 
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which gives Ii{2z) / Iq{2z) —* 1 and therefore 

(iV±)c.e. = (iV±),.c 



Z . 



(8) 



(9) 
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Using the series expansion one gets [7j for small systems {z <^ 1): 



In{2z) = — + 



yn+2 



n\ (n + 1)! 
and consequently Ii{2z) / Io{2z) = z which results in 

(iv±)c.e. = z^ (iv±) 



+ O {z^+^) 



g.c.e. 



Z . 



(10) 
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The asymptotics of the mean multiplicity discussed above are clearly seen in Fig. [H 

4. Scaled variance. An useful measure of fluctuations of any variable X is the ratio 
of its variance V{X) = (X^) — (X)^ to its mean value {X), referred here as the scaled 
variance: 



{X') - {xy 



;i2) 



Note, that = 1 for the Poisson distribution. Thus, to study the fluctuations of charged 
particles the second moment of the multiplicity distribution (N^) has to be calculated. 
In the g.c.e. (jT)) and c.e. © one finds: 



{Nl)c.e. 



Jg.c.e 
1 



d 



d 



g.c.e. 



d Z,, 



A±=l 



A±=l 



Z + z- 

h{2z) 
h{2z) 



+ z' 



h{2z) 
h{2z) 



(13) 
(14) 



The corresponding scaled variances are: 



UJ 



± 

g.c.e. 



{Nl)g.c.e. - {N^? 



g.c.e. 



{N±)g.c.e. 

{Niu, - (iV±)L. 



h{2z) h{2z) 



Io{2z) h{2z)\ 



(15) 
(16) 



Using Eqs. (jH)) and (fTUI) the asymptotic behaviour of tu^g for both z — > and z oo can 
be found. The c.e. fluctuations measured in terms of u are equal to those in the g.c.e. for 
the small system (2; ^ 1) (another variable to treat the fluctuations in the small systems 
is discussed in Appendix): 



UJ. 
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± 

g.c.e 



(17) 



For large systems {z ^ 1) the scaled variance for the c.e. is two times smaller than the 
scaled variance for the g.c.e. : 
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g.c.e. 



The dependence of the scaled variance calculated within the c.e and g.c.e. on z is shown 
in Fis;. 2. 
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Figure 2: The scaled variances of N± calculated within the g.c.e. , w^^.e. 
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The scaled variance shows a very different behavior than the mean multiplicity. In the 
limit of small z the ratio of the results for c.e. and g.c.e. approaches zero for the mean 
multiplicity (Fig. 1) and one for the scaled variance (Fig. 2). On the other hand in the 
large z limit the mean multiplicity ratio approaches one and the scaled variance ratio 
0.5. Thus in the case of fluctuations the canonical and grand canonical ensembles are not 
equivalent. 

5. Multiplicity distribution. In the g.c.e. the multiplicity distribution of (and 
A^_) is equal to the Poisson one: 



Ar_=0 



exTp[—z) ■ 



Z„ce ^ A^+! A^-! 



(19) 



whereas the corresponding distribution in the c.e. Q is: 



Af_=0 
1 



Zee ^ NJ. NJ 



5 {N+ - N_) 



(20) 



Io{2z) \N^\ 

As an example, the distributions in g.c.e. and c.e. are plotted in Figs. 3 and 4 for z = 0.5 
(the small system) and z = 10 (the large system), respectively. 
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Figure 3: Multiplicity distributions Pc.e.(A^±) ^ and Pg.c.e.iN±) ^ for z = 0.5. 
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Figure 4: Multiplicity distributions Pc.e.(iV±) (EOl) and Pg.c.e.{N±) ^ ior z = 10. 
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As expected from the previous discussion, the c.e. distribution (j^Uj) is narrower (the 
variance is smaller) than the g.c.e. one (fTU|) . This result is valid for both the large 
(2; ^ 1) and the small (2; -C 1) system. On the other hand, the average value of N± is 
smaller in the c.e. than in the g.c.e. for small z. It results in cj^g —>■ ujf^f, = 1 at 
z — > 0. Moreover, for {N±) <C 1 one can easily demonstrate that u;^ = 1 for any P{N±) 
distribution if the conditions -P(O) ^ -P(l) ^ -P(^) (with k > 2) are satisfied. Indeed, in 
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this limit one can neglect all P{N±) for N± > 2 which results in: 

, _ (iV|) - {N^)' ^ P(l) ■ 1' - [P(l) ■ 1]^ ^ 

as -P(l) = {N±) ^ 1. In the large volume limit, see Fig. 4, the mean values of the c.e. 
and g.c.e. distributions become equal, but the c.e. distribution is narrower than the g.c.e. 
one. 

6. Total multiplicity of charged particles. The total multiplicity of charged 
particles is defined as N^h = + N_. Its average in the g.c.e. and c.e. reads: 

{N,h)g.c.e. = {N+ + N_ = (iV+),.e.e. + {N.)g,,,e. = 2z , (22) 

{N,h)c.e. = ( iV+ + iV_ ),,. = (iV+),.e. + {N^)c.e. = 2z ^-^^ . (23) 

Io{2z} 

In the g.c.e. one finds: 

{N^h)9.c.e. = (iV2 +2iV+iV_ + iV!),.,.,. = {Nl)g,,,,, + 2{N+)g,,,,,{N.)g.,.e.+ (24) 
+ {Nl)g.c.e. = + z + 2z^ + z^ + z = 4z^ + 2z , 



and consequently the scaled variance of Nch in the g.c.e. is: 



(iVi)p.c.e. - {Nch)1,,,, ^ Az^ + 2z - (2zf 



The result ()25|1 also follows from explicit expression on the probability distribution of Nch 
in the g.c.e. : 



oo oo 



P,.c.e.(iVch) = 5^ 5^ P9.c.e. iV-) " 5 [Nch - (iV+ + iV-)] (26) 

N+ N-=0 



^oce ^ ^ A^+! A^-! iVc/»! 

Thus distributions of Nch and A^-t are Poissonian in the g.c.e. . In the c.e. the negatively 
and positively charged particles are correlated, (A^+ ■ N_)c.e. 7^ {N+)c.e. ■ (A^-)c.e.- The 
correlation term reads: 

Using Eqs. (fT^ and P7jl one obtains the scaled variance of Nch in the c.e. : 



.ch _ (A^i)c.e. - (iVch)' 



c.e. 



1 + 



l2{2z) + Ioi2z) _ ^ h{2z) 
h{2z) Ioi2z)_ 



i2t 



The scaled variances u':^^ ^ and u^^^ as functions of z are shown in Fig. together with 
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Figure 5: The scaled variances uj\ 
functions of z. 
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From Eqs. (|TT)|l and and the recurrence relation /o(2z) = hi^z) + Ii{2z)/z 7\ it 
follows that tu^ g = 2co'^g , i.e. the relative variance of total charge multiplicity Nch is two 
times larger than the one of N±. This is because Net = 2N^ = 2A^„ in each microscopic 
state allowed by an exact charge conservation. One obtains a similar result for the case of 
particle production via decay of neutral resonances, e.g., p° — ^ tt"*" + 7r~. The distributions 
of TT^ and 7c~ coincide with the p° distribution, and consequently uj^ = uj, where u is the 
scaled variance of the distribution of p°. But because N^h = 2Np one gets lo'^^ = 2a;. 

Probability distribution of N^h in the c.e. reads: 



Pc.e.{N,h) = J2Y1 (^+' ■ ^ [^-^ - + ^-)] 

N+ N-=0 

1 



(29) 



h{2z) 
1 

hi2z) 



N+=0 N-=0 

(iVch/2)! 



5(iV+-iV_).5[iV,,-(iV+ + iV_)] 



It coincides, of course, with Pc.e.(^+) (1201) at = Nch/2. As an example , the probability 
distributions Pg.c.e.{Nch) PHjl and Pee. (Nch) (1201) are shown for 2; = 0.5 (the small system) 
and for z = 10 (the large system) in Figs. 6 and 7, respectively. Only even multiplicities 
Nch = 0, 2, 4 . . . are allowed in the c.e. because of an exact charge conservation. For 
the small system {z < 1) the uo"'' reads (both Pg.c.e.{Nch = 1) < 1 and Pc.e.iNch = 2) -C 1 
at z < 1): 

t2 



ch 

g.c.e. 



UJ. 



ch 



P,.,.e.(l) ■ I' - [Pg.c.e.(l) • 1]^ ^ 

^..c.e.(l) -1 
P,.e.(2)-2^ - [Pc.e.{2)-2f - [Pc.e.{2)-2f 
Pc.e.i'2)-2 



(30) 
(31) 
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Figure 6: Multiplicity distributions of N^h for z = 0.5 in the g.c.e. and c.e. 
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Figure 7: Multiplicity distributions of Nch for 2; = 10 in the g.c.e. and c.e. . 
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In the large z limit the average number of charge particles {Net) and its scaled variance 
in the g.c.e. , Eqs. and (123), are equal to those in the c.e. , Eqs. and 
. Nevertheless the corresponding probability distributions are different, see Fig. 7. 
This is because all odd multiplicities are excluded in c.e. as a consequence of the charge 
conservation. The relation between Pg,c.e.{^ch) and Pc.e.{Nch) for the large 

system {z ^ 1) can be established as follows. Let us introduce the probability distribution 
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P*{Nch) defined as 



, 



0, 2, 4, 

1, 3, 5, 



where the constant C is given by a normahzation condition 

1 = J2 P*{Nch) = C- Yl Pg-ceiNch) 

Af,;, =0,2,4,... Ar,h =0,2,4,... 

(2z)2" 



(32) 
(33) 



(34) 



C-exp(-22) ^ 



n=0 



(2n)! 



C ■ exp(— 2^;) cosh(2z) 



Using Eq. (|34j) one gets C = 2 • [1 + exp(— 4^;)] ^ = 2 for z ^ 1. The origin of the result 
C = 2 is the fact that 



Pg.c.e.{N^h + l) = Pg.cANc 



2z 



ch 



iVch + 1 



Pg.ceiN. 



ch 



(35) 



for Nch close to its average value {Nch) 



g.c.e. 



2z ^ 1, i.e. if the odd numbers N, 



ch 



1, 3, 5, . . . are forbidden the probabilities Pg.c.e.i^ch) for the even numbers Net = 0, 2, 4, . . . 
should be approximately doubled to have a correct normalization for P*{Nch) 
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Figure 8: Multiplicity distributions PceX^ch) ^ and P*iNch) ^ for z = 10. 

Using the Stirling formula, n! = ri"'e~"V27m, valid for n S> 1, one finds that Pc.e.{Nch) — 
P*{Nch) for Ai'c/i close to its average value equal to 2^; ^ 1. Both distributions are plotted 
in Fig. 8 for a comparison. 

7. Limited kinematical acceptance. In the experimental study of nuclear collisions 
at high energies only a fraction of all produced particles which carry conserved charges 
is registered. Thus the multiplicity distribution of the measured particles is expected to 
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be different from the distribution of all produced particles. Within c.e. the effect of the 
limited kinematical acceptance (the acceptance in the momentum space) can be taken into 
account introducing a probabihty q that a single particle is registered. Because in c.e. 
particles are uncorrelated in momentum space the multiplicity distribution of accepted 
particles for a fixed number of produced particles N is given by the binomial distribution: 

Paccin, N) = - qf~- ■ ,,,, ' . (36) 

nl[i\ — n)\ 

Consequently one gets: 

n = qN, ^ -n^ = g(l -g) iV , (37) 

where {k = 1, 2) 



N 



Y.n" Pacc{n,N) . (38) 



n=0 



Introducing the probability distribution P{N) the first two moments of the distribution 
of accepted particles can be calculated: 



oo N 



{n)acc = $^i^(iV) $^nP,,,(n+,iV) = q-{N), (39) 

oo N 

{n')acc = $^P(iV) 5^n2p_(n,iV) = q'.{N') + q{l - q) ■ (N) , (40) 

N=0 n=0 

where {k = 1, 2) 

oo 

(A^'^) = Pi^) ■ (41) 

N=0 

Finally, the scaled variance for the accepted particles can be obtained: 

Uacc = q-^ + (1 - g) , (42) 

where uj in Eq. (j42j) is the scaled variance of the P{N) distribution. Assuming that P{N) 
corresponds to the c.e. , one finds from Eq. (j42p the scaled variance for the accepted 
particles in the c.e. tu^^^ = 1 for g <^ 1 and u^^^ = cu^g for q = 1. 

These limiting behaviour agrees with the expectations. In the large acceptance limit 
(g ~ 1) the distribution of measured particles approaches the distribution in the full 
acceptance. For a very small acceptance (g ~ 0) the measured distribution approaches 
the Poisson one independent of the shape of the distribution in the full acceptance. 

8. Summary. The particle number fluctuations have been considered within canonical 
ensemble for a system with zero net charge. The results are compared to those in the 
grand canonical ensemble where only the mean value of charge is required to be zero. In 
the large volume limit the fluctuations in c.e. are found to be different from those in the 
^f.c.e. . Thus the well known equivalence of both ensembles for the mean quantities is 
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not valid for the fluctuations. The scaled variance of the multiplicity distribution of same 
charge particles is calculated to be 0.5 in the c.e. and it is two times smaller than the 
scaled variance in g.c.e . These results may be relevant for the analysis of fluctuations 
in high energy nuclear collisions. In view of this the influence of the limited kinematical 
acceptance on multiplicity fluctuations also have been discussed. 

In this work the influence of the electric charge conservation was discussed. However, 
other material conservation laws, e.g. baryon number, strangeness or charm, can be 
treated within the same scheme. An extension of this work for a non-zero value of the 
conserved charge and several species of charged particles as well as an influence of an exact 
charge conservation on the energy fluctuations in the c.e. will be presented elsewhere. 
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Appendix. A variable = {X{X — 1))/ {X)"^ , was used for study of the fluctuations 
of N± and Nch in small systems {z ^ 1) 8J. In the g.c.e. (i.e. for the Poisson distribution 
of N± and Nch) one gets (-^2*^) = (-^2^'^) ce ~ whereas in the c.e. one finds (see 
also Fig. ©): 

.p±x ^ (iV± (N^ - _ Ioi2z).hi2z) ^ 1 ^ 

^ ' - (Ar±)L. " l!{2z) 2 + 6 ' 

^ (iVch {N,h - l))c.e. ^ Il{2z) + h{2z) ■ h{2z) ^ l_ 
^ ' ^'^■^ " {N,h)le. 2If{2z) 2z^ ■ 




Figure 9: The fluctuation measure F2 as a function of z in the c.e. 
solid lines indicate (-^2"^)^^ (03) and {F2^)^^ (EH), respectively. 
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In the large volume limit F2 1 and Fg'^'^ 1 because (A^±)^ 3> (A^±) ^ 1 (the 
same for Net)- Therefore, this measure is not suitable for a study of the particle number 
fluctuations in the large systems. 
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